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Abstract 

We present a theory of compatible differential constraints of a hydrodynamic hierarchy 
of infinite-dimensional systems, ft provides a convenient point of view for studying 
and formulating integrability properties and it reveals some hidden structures of the 
theory of intcgrable systems. Illustrative examples and new integrable models are 
exhibited. 

1 Introduction 

In this paper we are concerned with the hierarchy of evolution equations 
dG 

— = (Ai,G), G = G(X,x,t), t:=(h,t 2 ,...), i > 1, (1.1) 
where (U, V) := UV X — U X V, and the function G is assumed to posses an expansion 

G = go{x,t) H 1 ^ 1 ) A ^ oo, (1.2) 

We consider two different forms of (|l.lj) corresponding to the cases ffo = 1 and go ^ const, 
which will be henceforth referred to as the normalized G?W and the Schwarzian G^ 
hierarchies, respectively. They are determined by two different definitions of the functions 
Ai = Ai{X,x,t); namely, 

A (n) = V + • • • + gi ^X + 9i , 4 s) = g X l + ■■■ + 9l ^X, 
respectively. In both cases we may write 
Ai = (A*G9+, i > 1, 
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where the projectors P : A — > A + acting on power series A = ^ a n X n is defined as 

n 

A (n) := p (n) (A) = ^ anX n i j£) . = p« (A) = ^a n \ n . 

n>0 n>\ 



In terms of the coefficients of the expansion of G, the hierarchy becomes a set 

infinite-dimensional systems of hydrodynamic type 

n 

dign = S ^{gn-k,gi+k), n > 0, (normalized case), 
k=i 

n 

di9n = S ^(gn-k,9i+k), n > 0, (Schwarzian case). 

k=0 

We will show below that the case g$ ^ const can be considered as the Schwarzian form of 
the case go = 1 and that a transformation connecting both cases exists. 

Although equations similar to were considered in the literature ^|2j|2lE] (starting 
with pQ), the normalized case of the hierarchy was recently reconsidered in [3] from 
the point view of the reduction theory of integrable systems. It was in this context where 
the notion of differential constraint was introduced. One of the main aims of the present 
work is to show that this notion is a useful one in the theory of integrable systems. 

For several purposes it is convenient to use instead of G the function 

1 /ii ho 

ho — — , tl\ — ii2 o - n, V-L-oJ 

90 go 95 go 

It allows us to write in the alternative form 

which implies that the coefficients of the expansion of H in powers of A supply an infinite 
set of conservation laws for 

The function H is particularly useful to deal with the reductions of (jl.lj) determined 
by differential constraints. Let us illustrate this fact by considering the first flow t := t\ 
of the Schwarzian hierarchy, which according to (|1.4|) reads 

h n>t = h n+lx , h n :=^ } -, n> 0. (1.5) 

ho 

Let us look for a differential constraint of the form 

hi = B(h ,h 0tX ), (1.6) 
allowing to reduce the first equation in 1)1. 5|) to 



ho,t = D x B(h , ho, x ), 



(1.7) 
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where D x stands for the total differentiation with respect to x. The differential con- 
straint Q1.6B should be compatible with the whole hierarchy. In particular, the second 
equation in (|1,5|) implies 

A(M = D t (h B(h ,h 0iX )) G ImD x . (1.8) 

Here Dt is the total differentiation operator with respect to t and lmD x denotes the range 
of D x acting on the set of polynomials in the variables /io, ho )X , ho ;XX , .... Now one easily 
finds that 

D t (hi) = -(hoBnB^hl^ + ■■■ , 

where B\\ is the second derivative of B with respect to the argument ho tX . Thus (|1.8I) 
yields B\± = and therefore the constraint (|1.6|) should be of the form 

h = b(h ) + a(h )h QiX . 

One can directly prove in this way that the only constraint Q1.6JI compatible with (|1.5j0 
is 

(logh ) x = h 1 . (1.9) 
The corresponding reduced evolution equation Q1.7JI is given by 

V + ( j- ) = °- 

\ h 0/xx 

Several interesting generalizations of this type of constraints will be discussed in this 
paper. For example it will be proved that the constraints 

(log ho) x = h m , m > 1, (1.10) 

and 

{ho,x}= [H 2 ] n , n>l, (1.11) 
are compatible with (|1.5|) . Here 

{a ' x} ■= \ (t! 



1 &XX 

2~a~" 



and [-ff 2 ] n stands for the coefficient at X~ n of the expansion of H 2 in powers of A. 

The paper is organized as follows. In Section 2 we deal with the structure of symme- 
tries of the hierarchies (|l,lj) and reformulate some basic results. A particular attention 
is devoted to applying the classical Liouville method of solution for solving the polyno- 
mial reductions of (|1.1|) . Some apparently new results concerning r-functions like objects 
(Theorem 2) and a map -» (Th eorem 3) are also presented. Section 3 is con- 
cerned with the use of the method of differential constraints to characterize nonlinear 
integrable models from (jl.ljl . It is shown how the structure of wide classes of these models 
as well as simple strategies for finding solutions can be conveniently described by this 
method. Some new integrable models, including a 2 + 1-dimensional one (|3.22|) are exhibi- 
ted. 
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2 Symmetries 

We first reformulate some already known basic properties about the symmetries of the 
hierarchy (jl.ljl . Our starting point is the following theorem [S] 

Theorem 1. The evolutionary flows defined by form a commuting family. 

Proof. By using the Jacobi identity for the wronskian operation (U, V) we have 

&G___&G_ = /d_A L _d_A i 

dudtj dtjdu \dti dtj Kvih 

^,From Q1.1J1 and taking into account that Ai = (A J G)+ it follows 

^-^ = ((A i ,X j G)-(A j ,X i G)) + . (2.1) 

On the other hand 

{Aj,X l G) = -(X l G,X j G-(X j G)-) = (X i G, (A J G)_), 
where (X j G)- := X j G - (X j G)+. Hence 
dAi dA 



ti 

dti dtj 



{(A h XiG) - {X i G, (A'G)_» + = {(A,X^G) - (A l ,(X^G)^))_ 
(Ai,Aj}+ = (Ai,Aj) ; 



which proves the statement. ■ 

The hydrodynamic character of the equations of the hierarchy has important 

consequences. Thus, it follows that if A = Ao is either a zero or a pole of G then 

— = MX ,x,t) — . 

It means that zeros and poles of G are Riemann invariants of Ijl.ljl . This is a useful prop- 
erty because an obvious reduction of arises when only a finite number of coefficients 
in the expansion 1)1. 2|) of G are different from zero 



G^ = l + | + ... + fJ, G (s)= go + !l + ... + i» r ±. (2.2) 

In this case it is possible to perform the integration of in terms of Riemann invariants. 
Let us illustrate this feature for the normalized case. We can rewrite G in the form 

1 N 

G= w U(X + ll (x,t)), (2.3) 

i=l 

where ji(x,t) are the zeros of G. Thus, under the change of dependent variables 
9 ■= {9i,---,9n) ^7 := (7l> • • • > 7iv), 
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the system (jl.lj) reduces to 

dnli = Slni{l)dxli, n > 1. (2.4) 

Here Q ni := A n \ x= - 1% are given by 

:-l) n E7i 1 7 J v7i„ if n = l,2,...,iV-l, 
fi™ = { (2.5) 
if n > iV. 



The equations (|2.4|) ~ (|2.5|) form a set of N — 1 weakly nonlinear hydrodynamic systems 
of Dubrovin type 2 which can be integrated by means of a version of the classical Liouville 
method [£] proposed by Tsarev (see also |Hj). It basically consists in finding a vector 
field D x acting on the new dependent variables 7 

D x n=Xi{% i = l,...,N, (2.6) 

in such a way that it commutes with the fields dj, j = 1, . . . , N — 1. This is verified 
provided 

dlnXt 1 

i^j. (2.7) 



d 7i 7 j - 7j 
By applying Tsarev's method (see |S] ) one finds 
Lemma 1. ^4 general solution of (|2.7|) 

n(7i -7j) 

where dj(£) are TV arbitrary functions. 

As a consequence we can reformulate ()2.4j) as a set of iV dynamical systems for 7$ 

3 n7 i = VLniXi, n = 1, . . . , N - 1, 

= Xi. (2.9) 

Equivalently, we have 

^ AT-l \ 

dx+ J2 n nidt n ) , i = l,...,N. (2.10) 



d7i 



«i(7i) ]l(7i-7j) v „ , , 
One can explicitly invert Q2.1U|) (see 0) and find 

N N-n-1 

(-l)"di n = ^^- T ^d 7 i, n = 0,...,iV-l, (2.11) 

i=l a n7«J 

where we are denoting to := x. Therefore, we conclude that the general solution of 1)2.41) 
is determined by the following system of N implicit relations 

(-l) n t n = ^2 Vf^' n = 0,...,N-l. (2.12) 

i=i a «l7«J 

The theorems which follow state some apparently new results on the structure of the 
hierarchy (|1.1|) 
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Theorem 2. The differential forms 

u\ = gi d x + gi dt\ H 1- g n d t n+ \ + ■ • • , {normalized case), 

1 

0J2 = — (da; — gi dt% H — g n dt n + ■ • • ) , (Schwarzian case), 

9o 

are closed. 

Proof. Let us consider the normalized case. Notice that 
= ^Res(A i G)dt i , 

i>0 

where Res(^4(A)) stands for the coefficient of A -1 of Laurent series ^4(A) in A. Moreover, 
from ifLTT) 



Res(A*G) = Res (^^) = Res((A,-, A;G}). (2.13) 



On the other hand 

Res ((A,-, XiG)) = Res({X j G - (AG)_, AjG)) = -Res(((A J G)_, (A,G)+)) 
= Res((A 4 , (AjG)_)) = Res((Ai, A.-G}), 

so that ()2.13|) implies 

4 Res(A ' G) = l Res(A,G) ' 

which proves that u>i is closed. 

The proof for the Schwarzian case can be done similarly by writing 



L02 

9o 



L_y Res (V^-W 

70 ^ V W 



and by taking into account that 

d 1 d oa , 

W- = -7f- ■ 2 ' 14 
dti 9o ox g 



As a consequence of this theorem there exist two functions qi = qi(x,t) verifying 

Ui = d qi , i = 1,2, (2.15) 

so that we can rewrite the hierarchy (jl.l|) as a system of partial differential equations 
for q { . 



Towards a Theory of Differential Constraints of a Hydrodynamic Hierarchy 



235 



Theorem 3. Let = G^ s \\,x,t) be a solution of the Schwarzian hierarchy (jl.lj) . // 
we define x = q2(x,t), where q2 satisfies (|2.15j) . then 

G<">(A,s,t) == G?W( r A '^ t) , (2.16) 

solves the normalized hierarchy (jl.lj) . 
Proof. Observe that 

dx 1 dx 1 • > j 

Hence 

ox Oti 
Therefore, by taking (j2.14j) into account, we have 



<9t; +9t 'dxJ g g 

where 

/t ( n ) \ i i i A \ i * - 9i 

A\ —AH hft-iA + ft, ft := — , 

go 

and (£/, F) f := Ud&V - Vd £ U. U 

Remark. This theorem shows that a correspondence G^ s ' — > G^ n ) between the Schwarzian 
and normalized hierarchies exists. Reciprocally, a transformation G^ — ► G^ can also be 
defined provided that the solution G^ = G^ (X,x,t) does not vanish at A = 0. To prove 
this statement let us define 

1 

u :- 



and the differential form 



u) := u dx + y^ft dtj , 



i>l 

where g n are the coefficients of the expansion of in powers of A. From (jl.lj) it follows 
at once that 

diU = (ugi) x . 

Moreover, from the identities 
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and taking into account that gi = ^4j|A=0) we have 

djQi - dig-j = (gj,gi), 
and therefore 

dj(ugi) = di(ugj), 

which proves that uj is a closed form. It is now straightforward to demonstrate that the 
function 

G {s) {X,x,t) := u(x,t)G {n) (X,x,t), 

where 

x := q(x, t), uj = dq, 
verifies the Schwarzian hierarchy (jl.lj) . 

3 Reductions and differential constraints 

In [3] wide classes of reductions of the normalized hierarchy were introduced by using 
the notion of compatible differential constraints. They can be easily generalized to the 
Schwarzian hierarchy as well. We will consider here reductions defined by second-order 
differential constraints 

2G XX G -G 2 X + 4a(A) - 4£/(A, x, t)G 2 = 0, (3.1) 

where a(A) is an arbitrary function and 

U(\,x,t):=(^f) . (3.2) 

Requiring compatibility between IJl.lj) and (|3.1j) yields 

d n U = ~—A njXXX + 2UA n x + U x A n , n > 1. (3-3) 

We notice that in terms of the generating function H for the conservation laws, equa- 
tions (|3.1|) and (|3.2|) can be written as 

{H,x} + a(X)H 2 = U(\,x,t), U=(a{\)H 2 ) + , (3.4) 

where we are using the Schwarzian derivation operation 

3 H 2 . 1 H xx 



{H,x} r- 



4H 2 2 H 



The equations determine integrable hierarchies associated with Schrodinger spectral 
problems with energy-dependent potentials. Thus, for o(A) = A n one finds the generalized 
KdV hierarchies for energy-dependent potentials of the form 



n-l 



U(X,x,t) := X n + X l Ui(x,t), (normalized case), (3-5) 



i=0 
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and 

n 

U(X,x,t) :=^^ j X i Uj(x,t), (Schwarzian case). (3-6) 



i=i 



For the normalized case G = G^ the simplest choices n = 1 and n = 2 lead to the KdV 
and the Zakharov-Shabat hierarchies, respectively. We observe that for the Schwarzian 
case the choice a(A) = A n in 1)3.4)) leads to the constraint p.llj) . 

Theorem 4. The differential constraint (|3.4j) is invariant under the transformation 
_> jj{n) determined by equation (j2,16|) . That is to say, if satisfies 

(ff (s) ,i} + a(A)ff (s)2 = (7 (s) (A,x,t), 

then 

x} + a(X)H^ 2 = U( n \\,x,t). 
The corresponding transformation law for the potential function is 

U^=g 2 U^ + {g (x),x}. (3.7) 
Proof. If we define := H^/h , then it follows at once that 

{H^\x}=hl^{H^,x} + \ h ^ lk °^ 

On the other hand 



hi 2 h 



■2 



1-1 1 ho,£x 1 ^0,£ 
x — 



ho' J 2 h 4 hi ' 
Thus, it is straightforward to get 



{ffM.x} + a(X)H^ 2 = U {n \ C/ (n) := ^W + jJ-^J. 



(3.8) 



A complete and convenient description of the general class of the energy- dependent 
hierarchies can be formulated by taking advantage of the identities (|3,4j) . The following 
examples illustrate the fact that, in addition to the standard hierarchies of Zakharov- 
Shabat and KdV types, many other integrable models may be analyzed within the theory 
of the second order differential constraints of the Schwarzian hierarchy 1)1. 1J) . 

Example 1. For a(A) = A and G = G^ s \ the constraint 1)3.4)) reads 

{H,x} + XH 2 = Am, 
and by identifying the coefficient of 1/ A we get 

{h ,x} + 2h hx = 0. 
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Hence, by taking into account that in terms of the function q := q%(x,t) of (|2.15|) 

h = — = q x , hi = = q x q t , (t := h), (3.9) 
we deduce at once the Schwarzian form of the KdV equation 

2 

2q t q x = {q x ,x} o -8 ft + 2^£ - 3% = 0, (3.10) 

Qx Qx 

which, written in terms of go = 1/ho, reduces to the Harry-Dym equation 
1 3 

9o,t — -r9o9o,xxx- (3.11) 

Example 2. For a(A) = A 2 and G = G^ s \ the identity (|3~1)) reads 

{H,x} + X 2 H 2 = X 2 u 2 + Xu 1 , 
and by identifying the constant term in the expansion in powers of A we get 

{h ,x} + h 2 + 2h h 2 = 0. 
Now, from (jl.4j) we have 

Hence, by taking into account (|3.9fl we find 

/ \ / 1 Qxxx 3 q xx 1 2 \ /, 10 \ 

This nonlinear equation admits a Lagrangian given by 

QtQx + l^Adxdt, (3.13) 
Qx J 

or, equivalently, under the point transformation q(u(y,t),t) = y 

£ = U^kdydt. (3.14) 
If we denote now y — ► x, u — > g, the following equivalent Lagrangian formulation arises 
q t +q ^dxdt= l[( {qt + q r ? -2 q -^\dxdt 



(Qt + Q: 



\2 



xx I 



Qx 



dxdt. 



Thus, the model admits a Hamiltonian density 



1 2 2 r r . (Qt + Q 



\2 



H = q t Lq t ~ L = p x q x + -p q x , p := L gt , L 



XX I 



qt *-> — fxMx i .f Hxi V -~ ^qti • 2 

* Qx 
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which leads to the Hamiltonian equations 

qt + q xx + 2pql = 0, -p t + Pxx = (2p 2 q x ) x . 

It should be noticed that the last system is one of "canonical" forms of the isotropic 
Landau-Lifshitz model (see 

In [H] first-order differential constraints defining reductions of were also intro- 

duced. They adopt the form 

G x -a(X) -U(X,x,t)G = 0, (3.15) 

where a (A) is an arbitrary function and 

U(X,x,t):=-(^j . (3.16) 

Requiring compatibility between (|1,1|) and (|3.15|) gives 

+ (UA n ) x , n>l, (3.17) 

or, equivalently, 

d n U = -d x (U(X n G).)+, n>l,. (3.18) 
In terms of the function H equations (|3.15|) and (|3.16|) read 

H x + a(X)H 2 + U{X,x,t)H = 0, U = -(a(X)H) . (3.19) 



We notice that for the Schwarzian case if we set a (A) = — X m in (|3.15j) the constraint (jl.l()|) 
follows. 

First-order differential constraints determine in turn further reductions of the second- 
order differential constraints. Indeed, it follows easily that 

Theorem 5. Given a solution G of (|3,15|) then it satisfies 

2G XX G - G\ + a(Xf - 4C/(A, x, t)G 2 = 0, 
where U — > U is implemented by the Miura transformation 

U:=\U X + \U 2 . 

Among the reductions defined by first-order differential constraints one finds the Bur- 
gers hierarchy as well as a class of hierarchies of energy- dependent type. 

Example 3. If we impose on G = G^ the constraint 

G x = (A + u)G - A, u = h\ = -gi, 
we get the Burgers hierarchy 

d n u = d x {(d x + u) n u) . 
In particular the flow corresponding to t := t\ is the Burgers equation 

ut = u xx + 2uu x . 
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Example 4. For a constraint on G = of the form 

G x = (A 2 + Xui + uq)G - A 2 , Ul = h 1 = -g 1 , u = h 2 = gf - g 2 , (3.20) 
we get the following hierarchy for the energy- dependent potential U := A 2 + Xui + uq 
ui \ ( -d x \ nn+1 ( 1 



9n = n o ,o x n , (3-21) 



where 

R : = 

The flow corresponding to t := ti reads 



-ui d x + u 
1 



uo,t = ui,xx - (uoui) x , v>i,t = ~ {uf + u o) x , (3.22) 
which in terms of the function q := q\(x,t) of (|2.15j) reduces to 

Qtt ~ Qxxx ~ SqxQxt ~ QtQxx + (q x )l = 0. 

It is possible to characterize multidimensional integrable models by imposing compat- 
ible differential constraints to Ql.lj) . To illustrate this feature let us consider the con- 
straint ljM.2()|) for G = G^ and denote y := t\, t := t 2 . Thus we have the system of 
equations 

G x = (A 2 - \ gx + g\ - g 2 ) G - A 2 , (3.23a) 
G y = (X + gi ,G), (3.23b) 
G t = {X 2 + X 9l +g 2 ,G). (3.23c) 

Let us introduce the function 

fi h 1 
F :=uG = f + ^- + £ + ■■■ , u:~- 



X X 2 ' • G| A=0 
From (|3.23|) we get 

u x = (g 2 - gf) u, 

u y = (gi, x + 5i52 - gf) u , 

u t = {g2,x + g 2 - 92gl) u, 

so that the equations (|3.23f) written in terms of F read 

F x = (A 2 - X 9l ) F - A 2 / , (3.24a) 
F y = (A 3 - Xg 2 ) F - A 3 / - A/i, (3.24b) 
F t = (A 4 - A<? 3 ) F - A 4 /o - A 3 /i - A 2 / 2 . (3.24c) 

Now it is easy to deduce from Q3.24JI that 

F t = F xx + 2 9l F y - glF x , (3.25) 
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and, as a consequence, one finds that v := g\ and w = In it satisfy the nonlinear system of 
equations 

vt = v xx + 2vv y - v 2 v x + 2v x w x , 

wt = w xx + 2vw y — v 2 w x + w 2 . (3.26) 

We observe that equations (|3.24j) mean that equation (|3.25j) and, consequently, equa- 
tion (|3.2H|) admit separation of variables. If we expand F in powers of A then equa- 
tions (|3.24|) become 

fn,x = fn+2 ~ gifn+1, f-1 =0, f = U, (3.27a) 

fn,y = fn+3 ~ fl^/n+l) H = u 92, (3.27b) 

fn,t = fn+4 ~ 93fn+i, /~3 = ug 3 . (3.27c) 

By using (|3.27a|) all the coefficients f n can be expressed as differential polynomials in u 
and v = g\\ 

fo = u, fi = uv, f 2 = u x + uv, f 3 = (uv) x + vu x + uv 3 , 
Furthermore, the system (|3.27|) implies that polynomial reductions of the form 
fn = 0, V n > N, 

are Liouville integrable. Indeed, by noticing that Jn-i = const, we require N — 1 com- 
muting vector fields to integrate such a reduction. For example, for N = A the dynamical 
variables are /o, fi, fi and (j3.27j) can be written as 

c2 ■ 



d/o = (yf2 ~ dx + - ^J^j d y- a J^ dt > 



, f ( fih\ , /2 i h , . 
d/i = a 7— dx - — dy-a— dt, 

\ Jo / Jo Jo 

dj 2 = -a—dx- a— dy - — dt, 
Jo Jo Jo 

where a 3 := = const. It is now straightforward to eliminate dx, dy, dt and find 

V a 2a 2 J 
i f fo fih . /| 

dy = d 2- + T^ 

a a z 3cr 



,,_ , / /0/2 _ /1 /1/2 _ /2 "\ 
V « 2 2a 2 + a 3 4aV ' 

In this way the solution is determined by the following system 

fi fi 

fo = ay + xf 2 - — , fi = ax + —, 
2a 2a 

- 7 -fi + \a 2 xf 2 + a 3 / 2 + a 4 Qx 2 + 1) = 0. 

This example shows the complicated algebraic singularities exhibited by the solutions 
of (I3~2lfl) . 
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